Spin accumulation in ballistic Rashba bar 
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We propose an analytic model to study intrinsic spin polarization effect in a ballistic Rashba bar 
with two semi-infinite leads. The wave functions expanded with plane waves in Rashba bar are 
required to satisfy boundary conditions at both longitudinal and transverse interfaces. We find out- 
of-plane spin Hall accumulation effect can be induced in the Rashba bar even with large dimensions 
by injecting unpolarized current from the lead. The longitudinal in-plane spin Hall effect, however, 
becomes obscure in large-size sample. An interesting direction-flipping of the out-of-plane spin 
accumulation is predicted by altering the Rashba coupling strength. 

PACS Number: 72.25.Dc,73.23.Ad, 85.75.Nn 
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Introduction: An external electric field can be expected 
to induce a transverse spin current in non-magnetic semi- 
conductors or result in out-of-plane spin accumulation 
near the edges of semiconductor films due to the spin- 
orbit coupling (SOC) of electron. It is the so-called 
spin Hall effect (SHE) 01 . Very recently, two indepen- 
dent groups have reported the experimental observation 
of SHE[2j, y], giving evidence of the existence of sig- 
nificant SOC in semiconductors. Such new effect of 
SOC creates a new way to manipulate electron spins 
by means of an electric field, not magnetic field, that 
may have great potential applications in future spin- 
tronics. SHE was first proposed by D'yakonov and V. 
I. PereP in 19710, then HirschQ and Zhangj^. The 
existence of SHE effect in p-type semiconductors and 
n-type semiconductors in two-dimensional heterostruc- 
tures was predicted by Murakami et aZ.@] and Sinova et 
aLjjj, respectively. However, the microscopic origin of 
SHE still has controversial explanations. The discussion 
of extrinsic or intrinsic origin currently attracts many 
scientists. 0, 0, (HI E2, EE Q EE EE 113 This letter 
focuses on the study of spin accumulation in Rashba SOC 
system in ballistic region, since ballistic samples usually 
with small dimensions are always desirable for future 
electronic device applications. References [lil H^. Efl Elf 
investigated spin accumulation in ballistic system with 
Rashba interaction. Based on discrete tight-binding 
model and Landauer-Biittiker formula, Nikolic et a/. [181] 
numerically demonstrated opposite spin accumulation 
near the two edges of the two dimensional semiconduc- 
tor structure. The phenomena is qualitatively similar to 
the experimental observations^!!!- Very recently, Wang 
et al. PJfl reported current induced local spin polarization 
due to Rashba SOC in a two dimensional narrow strip, 
where no lead is considered to attach to the Rashba strip. 
They concluded that out-of-plane spin polarization is im- 
portant mainly in systems with mesoscopic sizes, and it 
appears not to be associated with the SHE in bulk sam- 
ples. 

In this paper, we propose an analytic model to study 
intrinsic spin polarization effect in a Rashba bar con- 



nected to two semi-infinite leads. The geometry studied 
is similar to that proposed by Nikolic et al. 18]. We, 
however, don't use the discrete tight-binding model, but 
a continuous model to find out its exact solution to 
Schrodinger equation. The method we developed can 
properly deal with the coupling between channels in the 
Rashba bar. Therefore, the behavior of spin accumula- 
tion as a function of the width of the Rashba bar can be 
studied systematically. We find spin accumulation can 
be induced in Rashba bar by injecting unpolarized cur- 
rent from the lead. With the increase of the width of 
the Rashba bar, the peaks of spin accumulation shift to- 
wards the two lateral edges and the oscillations within the 
bar decrease, indicating the polarization becomes more 
prominent in the ballistic bar. This conclusion is very 
different from that obtained by Wang et al.\v^. Mean- 
while, we also report an interesting direction-flipping of 
the out-of-plane spin accumulation by altering the spin- 
orbital coupling strength, which may be considered as a 
signature of SHE in ballistic region. 

Theoretical Approach: We consider a 2DEG Rashba 
bar with two opposite interfaces connected to two semi- 
infinite ideal leads respectively along x direction. To sim- 
plify the calculation (not losing the essential physics) , the 
leads and the Rashba bar are assumed to have the same 
widths in y direction. And open boundary condition is 
applied at the two lateral edges of the strip. An elec- 
tron wave is injected from the right lead to the left one 
and crossing the middle Rashba bar. We solve separately 
Schrodinger equations of the three regions. The Hamil- 
tonian of the structure is described by 
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where m* is effective mass of electrons and a the Rashba 
SOC strength that can be tuned by gate voltage[22j. 
a = in leads, and a / in the Rashba bar. 
V{y) is the confined potential in y direction. For 
the open boundary, it is simply taken as V(y) = 0, 
(0 < y < b);V(y) = oo, (y < 0;y > b), where 
b is the width of the strip. The total wave func- 
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tion contains two components, ip(x,y) — ipi(x,y)Q ) ) + 
if>2(x,y)(i), where ipi,ip2 can generally expressed 
as i>i(x,y) = J2k n ,k x c n( k x)Xn(y)e ik * x ,Tp2{x,y) = 
J2k n ^ D n( k x)Xn(y)e lk=rX , where Xn(y) = sm(k n y), n = 
1,2,3.... In the formula, k n — rm/b and n is the chan- 
nel number. They are the wave functions that sat- 
isfy the hard-wall boundary conditions in y direction, 
i))(x, y)\ y= Q. y= i) — 0. With above wave functions, the 
Schrodinger equation is simplified as: 

E n (k x )C n (k x ) + iu ^y ^ [S nn 'k x — /\ n ^ n /k n r }D n ,{k x ) = 0, (2) 
n' 

E n (k x )D n (k x ) - ia^2{5 nn >k x + A n , n <k n >]C n >(k x ) = 0, (3) 

n> 

where E n (k x )=e n (k x )-E, e n (k x )=h 2 (kl+kl)/2m*, 
n,n'=l,2,3...,N, and A„, n /=§ J* sin(fc n y) cos(k n >y)dy, 
(n =/= n'). N is the total number of channels existing in 
the Rashba bar and dependent on the width b of the 
strip and the energy E of incident electrons. At fixed 
energy E, the total number TV is large with large width 
b. The terms containing A n/rl / in Eq.J2J) and JjJJ indicate 
explicitly the mixing interactions between different sub- 
bands due to the Rashba coupling and the boundaries 
existing in y direction. Our investigation reveals that no 
spin accumulation can be observed for the case of single 
channel in the bar, while the accumulation appears 
in multi-channel cases. This implies the mixing terms 
between subbands play an essential role in forming the 
spin accumulation effect. 
From Eq.© and ®, 

we have a polynomial equation 
of k x with the highest power of AN. Therefore, AN solu- 
tions of k x ]) ( j = 1 -> AN) can be obtained from the 
equation for any fixed E and N. In the case of sin- 
gle channel, where the SOC lifts the degenerate spin up 
and down states at ±k x into 4 different wave vectors 

ik x i and ik X 2- In general, N subbands would produce 

(?) (i) 

AN different k x . Each k x corresponds to an eigenfunc- 
tion <j>^ which has a length of 2N. All the coefficients 
{C n (k x ^), D n (k x ^)} in $^ are determined up to a nor- 
malized constant. And the wave function in Rashba re- 
gion with eigenvalue E must be a linear combination of 
AN eigenf unctions $ w , j = 1, 2, • ■ • , AN. The combined 
coefficients can be denoted by coefficients A(k x ^). There- 
fore, in Rashba region, there are totally AN coefficients 
{A(k x ^)} to be determined. If the Rashba bar is at an 
equilibrium state, AN eigenfunctions in Rashba bar have 
equivalent weight. There will be no spin polarization 
in the sample since Rashba interaction does not change 
the time-reversal symmetry. The time-reverse symmetry, 
however, is broken in the bar if there is steady current 
flowing along x direction. And spin polarization will oc- 
cur in the sample. Suppose an unpolarized incident elec- 
tron wave injected from the right lead. There are 2N co- 
efficients describing reflected waves in the right lead with 



respect to N channels for spin up and down states. Simi- 
larly, additional 2N coefficients are introduced to express 
the transmitted waves in the left lead. Therefore, there 
are totally 8N coefficients to be determined. They can 
be solved by boundary conditions at the two interfaces. 

The interfaces between leads and Rashba bar de- 
noted as i = 0, a are simplified by a 5 potential|23| 
VqS(x — xq),xo ~ 0, a (a is the length of the Rashba 
bar). The continuous conditions of wave functions 
and particle current crossing each interface, say at 
x = 0, yield equations: tp L (x)\ x=0 - = -i/j mid (x)\ x=0 + 

a,nd^t^ L (x)\ x=0 - - f£^^mw(aOU=o+ = if^i(O), 

where Pl x = (h/i)d/dx, and P mld x = (h/i)d/dx - 
(am* /h)(j y |U It has been assumed that there's 

no spin-flip across the interfaces [2j|. By means of the 
orthonormality of Xn, each subband of spin-up or -down 
state would satisfy above boundary conditions. There- 
fore, for the structure with N channels, two interfaces 
along x direction will yield 8N equations, the same as 
the number of coefficients to be determined. All the co- 
efficients in the wave functions thus can be exactly solved. 
With the wave function obtained, the local spin polariza- 
tion (Si(x,y)) = ^(tp(x,y)\ai\il>{x,y)) is straightforward 
to be calculated in the Rashba region. 

Results and Discussion: In the calculation, we focus 
on the linear region and simply assume that the Fermi 
levels and the effective masses of electrons in Rashba bar 
and two leads are the same. They are taken as 10 meV 
and 0.04 m e , typical values in 2DEG 22, 25]. The cur- 
rent injection has been normalized to one non-polarized 
electron wave for each channel. Meanwhile, we do the 
statistical average for the spin polarization of each inci- 
dent electron wave with energy E in each channel. The 
spin polarizations of (S z ), (S x ) and (S y ) with different 
numbers of channel are demonstrated in Fig. I, where av- 
erage integrals along x direction are taken. Note that 
the scale in the horizontal axis is normalized in order to 
compare explicitly the trends. The channel numbers of 
3, 6, and 9 correspond to the widths of about 120, 210, 
and 300 nm, respectively. The length of Rashba bar in 
the direction of current flow is fixed at 210 nm. It can 
be found that these sizes are within the phase-coherent 
length (L$ < O.A^j, - I^)|H Eil. Ex p erimentally per- 
mitted value of 2.9xI0 -11 eVm|22 l |27 | is taken for the 
strength of Rashba coupling. In the case of 3 channels, 
(S z ) has two high peaks with inverse symmetry distri- 
bution along y direction and there are small oscillations 
in the middle of the sample. With the increase of the 
channel number, the two high peaks become narrow and 
shift gradually to the lateral edges with the inside oscil- 
lations tending to be smoothed out. The phenomenon 
of out-of-plane spin polarization with opposite sign near 
the two edges is very much similar to the prediction 
of spin Hall effect by HirschQ and recent experimen- 
tal observations 0|- We also calculated spin polariza- 
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FIG. 1: The average spin accumulation of (S z ) (a), (S x ) (b) 
and (S y ) (c) as a function of the lateral position in the Rashba 
bar. The spin polarizations varied with the number of channel 
are indicated with different curve styles. 



tion of the Rashba bar having 20 channels (about 630 
nm in y direction). The envelop of its accumulation has 
the same behavior and keeps very strong spin Hall ac- 
cumulation effect. Therefore, the out-of-plane spin accu- 
mulation becomes more prominent with the increase of 
the channel number in the sample, namely, the width of 
the sample. This tendency is different from the conclu- 
sion obtained by Wang et al. pjfl • They reported that the 
out-of-plane polarization increases at small width of sam- 
ple, decreases, however, at large width with L ~ 16/ kp 
(about 6 channels). From our result, no evidence is found 
to support that their such tendency[19j can also exist in 
the Rashba bar with two leads at least to the width of 
the bar with 20 channels. 

The in-plane spin accumulations of (S x ) and (S y ) are 
given in Fig. 1(b) and (c), respectively. The distribution 
of (S x ) as a function of transverse position shares the 
same feature of (S z ), also forming SHE effect 0, Q, El- 
With the increase of the channel number, the effect of 
(S x ) also becomes more obvious. For (S z ) and (S x ), the 
net total spin polarization is zero because of the inverse 
symmetry of the distribution along y direction. The spin 
accumulation of (S y ), however, has a uniform sign along 
the y direction, which gives a net total spin polarization. 
The net transverse in-plane magnetization induced from 
unpolarized current can be ascribed to the combined fac- 
tors of Rashba SOC and time-reversal symmetry break- 
ing in the structures studied. 

From Eq. (2) and (3), the behavior of spin polar- 
izations of (S~)< (S x ), and (S y ) with opposite Rashba 
strength (— a)[l8| can be obtained analytically. In Eq.(2) 
and (3), an opposite Rashba coefficient —a produces an 
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FIG. 2: Spin accumulation of (S z ) (a) and (S x ) (b) with 
two different Rashba strengths (2.9xl0 -12 eVm, dotted line; 
2.9xlO~ n eVm solid line). The interesting direction-flipping 
is observed for (S z ), while no flipping is obtained for (S x ). 



opposite set of D n (k J x ), while C n (k x ) is not varied. There- 
fore, V'll-a = V'lla and V^l-a = -ip2\a ( see E q.(2) and 
(3)), which directly yield (S z (x,y)}- a = (S z (x,y)) a and 
(S x ,y(x,y))- a = —(S x ,y(x,y)) a . When the direction of 
Rashba interaction is reversed, the in-plane spin accumu- 
lation flips the polarization directions, while the out-of- 
plane one does not. We demonstrate, however, that in 
a ballistic SOC strip, the sign of (S z ) can be flipped by 
modulating the SOC strength, as depicted in Fig. 2(a). 
The thick curve in Fig. 2 (a) shows the spin polarization 
of (S z ) with small Rashba strength. With the increase of 
the strength to a certain value, opposite peaks near the 
lateral edges emerge, and the amplitudes magnify grad- 
ually while the inside oscillations decrease. The flipping- 
value of a is revealed to be inversely proportional to the 
width of the sample. No flipping, however, is observed 
for (5a;) by modulating the value of Rashba strength in 
the whole experimental permitted range [i^ H3 (see Fig. 
2(b)). Since the Rashba interaction strength in a normal 
2DEG can be tuned through gate voltage, we suggest a 
direct experimental detection of the spin- flipping of (S z ) 
crossing a critical value of Rashba SOC. 

The spin polarizations for (S z (x, y)) and (S x (x,y)) 
with long length (about 660 nm) of the Rashba bar are 
demonstrated in Fig. 3 (a) and (b), respectively. With 3 
channels in the bar, very obvious period-like variation of 
the pattern along x is observed for both (S z ) and (S x ), 
giving rise to the spin precession in the sample due to 
the SOC. The rough period of the precession is found 
to be closely related to the characteristic Rashba length 
L so = The 

rough period will be shorter if there 
are more channels in the sample or larger SOC. In Fig. 
3(a), the sign of (S z ) keeps in period-like variation along 
x direction, while (5a;) flips its sign during the variation 
(see Fig. 3(b)). This featured pattern of (S x ), how- 
ever, is obscured with the increase of the width of the 
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FIG. 3: Contour plots of spin accumulation of out-of-plane 
(S z (x,y)) (a) and longitudinal in-plane (S x (x,y)) (b) with 
long length of rashba bar. There are 3 channels in the sample. 
The Rashba strength used is 2.9 xlO~ eVm. 



sample due to having more oscillations along x direction. 
Therefore, for large dimensions of ballistic bar, the out- 
of-plane SHE can survive, while the longitudinal in-plane 
spin Hall polarization may not be detectable in the case. 

The variation of spin-polarization pattern under op- 
posite bias voltage[2 is predicted. If the center of the 
Rashba bar is taken as the origin, the wave function 
should have following symmetry: {i)j(x,y)\Si\ip(x,y)) -v = 
{ip(x,y)\P t U t SiUP\ip(x,y)) v , where P is the rotation of 
7r angle around z axis in orbital space, and U is SU(2) 
rotation with same angle ir in spin space. It is straight- 
forward to know that U = ( t ^ ■ I • Therefore, following 



relations can be obtained for the case of reverse bias. 

{4>(x,y)\S z \ip(x,y))- V = (ip(-x,-y)\S z \4>(-x,-y)) v , 
(ip(x,y)\S x \tp(x,y))-v = -y)\S x \i>(-x, -y))v, 

{ip(x,y)\S y \ip{x,y))-v = -(^(-x,-y)\S y \^(-x,-y))v. 

Our numerical calculations have exactly confirmed these 
symmetric relations, which are not completely same as 
the ones reported by Nikolic et aZ. 18]. For (S z (x, y)) un- 
der reverse bias voltage, the pattern of spin polarization 
does not just flip its sign 2], but has inverse symmetry 
compared to the original one (see Fig. 3(a)). It is also 
an intrinsic feature of SHE in ballistic region. 

Summary: We studied intrinsic spin polarization ef- 
fect in a ballistic Rashba bar by proposing an analytic 
model. Both the out-of-plane and longitudinal in-plane 
spin Hall accumulation effects are observed. In the sam- 
ple with large dimensions, only the out-of-plane spin Hall 
accumulation can be observed obviously. The transverse 
in-plane spin accumulation shows same sign across the 
sample, giving net spin magnetization. This spin non- 
conservation can be ascribed to the Rashba interaction 
together with the time-reversal symmetry breaking in the 
structure. An interesting direction-flipping of the out-of- 



plane spin accumulation is found by altering the Rashba 
coupling strength. This phenomenon is suggested to be 
confirmed by experiment. Some symmetric relations of 
spin polarizations are predicted. 
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